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Summary. Various issues related to the direct detection of supersymmetric dark 
matter are reviewed. Such are: 1) Construction of supersymmetric models with a 
number of parameters, which are constrained from the data at low energies as well 
as cosmological observations. 2) A model for the nucleon, in particular the depen- 
dence on the nucleon cross section on quarks other than u and d. 3) A nuclear 
model, i.e. the nuclear form factor for the scalar interaction and the spin response 
function for the axial current. 4) Information about the density and the velocity 
distribution of the neutralino (halo model). Using the present experimental limits 
on the rates and proper inputs in 3)-4) we derive constraints in the nucleon cross 
section, which involves l)-2). Since the expected event rates are extremely low we 
consider some additional signatures of the neutralino nucleus interaction, such as 
the periodic behavior of the rates due to the motion of Earth (modulation effect), 
which, unfortunately, is characterized by a small amplitude. This leads us to ex- 
amine the possibility of suggesting directional experiments, which measure not only 
the energy of the recoiling nuclei but their direction as well. In these, albeit hard, 
experiments one can exploit two very characteristic signatures: a)large asymmetries 
and b) interesting modulation patterns. Furthermore we extended our study to in- 
clude evaluation of the rates for other than recoil searches such as: i) Transitions to 
excited states, ii) Detection of recoiling electrons produced during the neutralino- 
nucleus interaction and iii) Observation of hard X-rays following the de-excitation 
of the ionized atom. 



1 Introduction 

The combined MAXIMA-1 [T], BOOMERANG 0, DASI 0, COBE/DMR 
Cosmic Microwave Background (CMB) observations 0], the recent WMAP 
data [S] and SDSS 6 imply that the Universe is flat [7] and and that most of 
the matter in the Universe is dark, i.e. exotic. 

Q b = 0.044 ± 0.04, O m = 0.27 ± 0.04, Q A = 0.69 ± 0.08 
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for baryonic matter , cold dark matter and dark energy respectively. An analy- 
sis of a combination of SDSS and WMAP data yields j§] f2 m « 0.30±0.04(lcr). 
Crudely speaking and easy to remember 

Q b w 0.05, flcDM « 0.30, n A w 0.65 

Additional indirect information comes from the rotational curves [jj]. The ro- 
tational velocity of an object increases so long is surrounded but matter. Once 
outside matter the velocity of rotation drops as the square root of the dis- 
tance. Such observations are not possible in our own galaxy. The observations 
of other galaxies, similar to our own, indicate that the rotational velocities of 
objects outside the luminous matter do not drop. So there must be a halo of 
dark matter out there. Since the non exotic component cannot exceed 40% of 
the CDM |J], there is room for exotic WIMP's (Weakly Interacting Massive 
Particles). 

In fact the DAMA experiment |IQ] has claimed the observation of one signal 
in direct detection of a WIMP, which with better statistics has subsequently 
been interpreted as a modulation signal These data, however, if they are 
due to the coherent process, are not consistent with other recent experiments, 
see e.g. EDELWEISS and CDMS It could still be interpreted as due to 
the spin cross section, but with a new interpretation of the extracted nucleon 
cross section. The above developments are in line with particle physics consid- 
erations. Thus, in the currently favored supersymmetric (SUSY) extensions 
of the standard model, the most natural WIMP candidate is the LSP, i.e. 
the lightest supersymmetric particle. In the most favored scenarios the LSP 
can be simply described as a Majorana fermion, a linear combination of the 
neutral components of the gauginos and Higgsinos |%]-|13|. 

Since this particle is expected to be very massive, m x > 30GeV, and 
extremely non relativistic with average kinetic energy T < lOOKeV, it can 
be directly detected mainly via the recoiling of a nucleus (A,Z) in the elastic 
scattering process: 

x + ( A) z) ^ x + (Azy (i) 

(x denotes the LSP). In order to compute the event rate one needs the fol- 
lowing ingredients: 

1. An effective Lagrangian at the elementary particle (quark) level obtained 
in the framework of supersymmetry [Hj, |14| and |13j . One starts with 
representative input in the restricted SUSY parameter space as described 
in the literature, e.g. Ellis et al |TS], Bottino et al, Kane et al., Castano 
et al. and Arnowitt et al ^1] as well as elsewhere ^J-^?]. We will not, 
however, elaborate on how one gets the needed parameters from super- 
symmetry, since this topic will be covered by another lecture in this school 
by professor Lahanas. For the reader's convenience, however, we will give 
a description in sec. of the basic SUSY ingredients needed to calculate 
LSP-nucleus scattering cross section. Our own SUSY input parameters 
can be found elsewhere [T%| . 
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2. A procedure in going from the quark to the nucleon level, i.e. a quark 
model for the nucleon. The results depend crucially on the content of 
the nucleon in quarks other than u and d. This is particularly true for 
the scalar couplings as well as the isoscalar axial coupling JH|-|2ni- Such 
topics will be discussed in sec. 0] 

3. computation of the relevant nuclear matrix elements j23"E21 using as 
reliable as possible many body nuclear wave functions. By putting as 
accurate nuclear physics input as possible, one will be able to constrain 
the SUSY parameters as much as possible. The situation is a bit simpler 
in the case of the scalar coupling, in which case one only needs the nuclear 
form factor. 

4. Convolution with the LSP velocity Distribution. To this end we will con- 
sider here Maxwell-Boltzmann B, (MB) velocity distributions, with an 
upper velocity cut off put in by hand. Other distributions are possible, 
such as non symmetric ones, like those of Drukier [23] and Green [24] . 
or non isothermal ones, e.g. those arising from late in-fall of dark matter 
into our galaxy, like Sikivie's caustic rings |25| . In any event in a proper 
treatment the velocity distribution ought to be consistent with the dark 
matter density in the context of the Eddington theory . 

After this we will specialize our study in the case of the nucleus 127 1, which 
is one of the most popular targets JU| > EH • 

Since the expected rates are extremely low or even undetectable with 
present techniques, one would like to exploit the characteristic signatures pro- 
vided by the reaction. Such are: 

1. The modulation effect, i.e the dependence of the event rate on the velocity 
of the Earth 

2. The directional event rate, which depends on the velocity of the sun 
around the galaxy as well as the the velocity of the Earth, has recently 
begun to appear feasible by the planned experiments 

3. Detection of other than nuclear recoils, such as 

• Detection of 7 rays following nuclear de-excitation, whenever possible 

133,1311. 

• Detection of ionization electrons produced directly in the LSP-nucleus 
collisions 

• Observations of hard X-rays produced 34] , when the inner shell elec- 
tron holes produced as above are filled. 

In all calculations we will, of course, include an appropriate nuclear form 
factor and take into account the influence on the rates of the detector energy 
cut off. We will present our results a function of the LSP mass, m x , in a way 
which can be easily understood by the experimentalists. 
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2 The Nature of the LSP 

Before proceeding with the construction of the effective Lagrangian we will 
briefly discuss the nature of the lightest supersymmetric particle (LSP) focus- 
ing on those ingredients which are of interest to dark matter. 

In currently favorable supergravity models the LSP is a linear combina- 
tion [H] of the neutral four fermions B,W3,Hi and H 2 which are the super- 
symmetric partners of the gauge bosons and and the Higgs scalars Hi 
and H 2 . Admixtures of s-neutrinos are expected to be negligible. 

In the above basis the mass-matrix takes the form 



/ Mi 

M 2 

-m z c/3S w m z Cf3C w 
V m z Sf3S w -m z C/3C w 



-m z C(iS w m z sps w ^ 
m z c f 3C w -m z Sf3C w 

-/! 

-H / 



(2) 



In the above expressions cw — cosOw, sw = sinBw, cp = cos/3, sp = sin/3, 
where tan(3 — (1)2) /{vi) is the ratio of the vacuum expectation values of the 
Higgs scalars H2 and Hi. /1 is a dimensionful coupling constant which is not 
specified by the theory (not even its sign). 

By diagonalizing the above matrix we obtain a set of eigenvalues rrij and 
the diagonalizing matrix Cy as follows 



/ B R \ 




/ XlB.\ 




( ^ L \ 




( XlL \ 


w 3R 

H\R 


= (Cg) 


X2R 
X3R 




HlL 




X2L 
X3L 


\H 2R ) 




\X4,R J 




\H 2 l) 




\X4lJ 



(3) 



with Cf- = C* e^ The phases are = 0,7r depending on the sign of the 
eigenmass. 

Another possibility to express the above results in photino-zino basis 7, Z 

via 

W3 = sinOwl — cosOwZ , Bq = cosOwl + sinOy/Z (4) 

In the absence of supersymmetry breaking (Mi = M 2 = M and /1 = 0) the 
photino is one of the eigenstates with mass M. One of the remaining eigen- 
states has a zero eigenvalue and is a linear combination of H\ and H 2 with 
mixing sin(3. In the presence of SUSY breaking terms the B, W3 basis is su- 
perior since the lowest eigenstate xi or LSP is primarily B. From our point of 
view the most important parameters are the mass m x of LSP and the mixing 
Cj\,j = 1, 2, 3, 4 which yield the xi content of the initial basis states. 
We are now in a position to find the interaction of xi with matter. We dis- 
tinguish three possibilities involving Z-exchange, s-quark exchange and Higgs 
exchange. 



On The Direct Detection of Dark Matter 5 

3 The Feynman Diagrams Entering the Direct Detection 
ofLSP. 

The diagrams involve Z-exchange, s-quark exchange and Higgs exchange. 
3.1 The Z-exchange contribution. 

This can arise from the interaction of Higgsinos with Z (see Fig.^) which can 
be read from Eq. C86 of Ref. [T3] 

L = —^-- A [H lia ^H 1R - Hil^Hil ~ (H2RJ»H2R - H 2L ^H 2L )]Z" (5) 
cosOw 4 

Using Eq. and the fact that for Majorana particles XJnX — 0, we obtain 




Fig. 1. The LSP-quark interaction mediated by Z-exchange. 



i=^^2(l^i| 2 -|C 4 i| 2 )xi7M75Xi^ (6) 
cosuw 4 



which leads to the effective 4-fermion interaction 

Leff = — ^— i2(|C 31 | 2 - |C 41 | 2 )(~ 9 2 1 2 X!7 M 75Xi)^ (7) 
cosOw 4 2cos6w q — 



where the extra factor of 2 comes from the Majorana nature of xi ■ The neutral 



hadronic current Jjf is given by 



J\ = -qi\{^sin 2 6 w - i(l - 75) - sin 2 9 w r 3 }q (8) 
at the nucleon level it can be written as 

Jf = -%a{ sm 2 9 w - g v (- - sin 2 6w)T 3 + ^9a7 5 t 3 }N (9) 
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Thus we can write 



L eff = -^(x i7 ^ 7 5 xi)Ja(Z) (10) 



where 



and 



UZ) = N lx [f*{Z) + ti(Z)r 3 + /°(Z) 75 + /A(Z) 75 r 3 ]7V (11) 



f v {Z) = 2(|C 31 | 2 - |C 4 i| 2 ) f Z 2 sm 2 6 w (12) 

m z — q z 

ti(Z) = -2(|C 31 | 2 - IQ1I 2 ) f Z 2 9v{\ ~ sinH w ) (13) 

m z — q z 2 

f A (Z) = , /A(Z) = 2(|C 31 | 2 -|C 41 | 2 )^^ \g A (14) 

m z — q z 2 

with gv — 1-0 and g A — 1.24. We can easily see that 

ti(z)/fv(z) = -gviTT^r - 1) ^ -i-i5 (15) 

Note that the suppression of this Z-exchange interaction compared to the 
ordinary neutral current interactions arises from the smallness of the mixing 
C 3 i and C41, a consequence of the fact that the Higgsinos are normally quite a 
bit heavier than the gauginos. Furthermore, the two Higgsinos tend to cancel 
each other. 

We should also mention that the vector contribution, the time component 
of which can lead to coherence, contributes only to order v/c~ 10 -3 due to 
the Majorana nature of the LSP. Thus to leading order only the axial current 
can contribute to the direct detection of the ncutralino. 



3.2 The s-quark Mediated Interaction 

The other interesting possibility arises from the other two components of %i, 
namely B and W3. Their corresponding couplings to s-quarks (see Fig- EJ) can 
be read from the appendix C4 of Ref. J3| They are 

L eff = -gV2{q L [T 3 W 3R - tan8 w (T 3 - Q)B R ]q L 

- tan9 w q R QB L q R } + HC (16) 

where q are the scalar quarks (SUSY partners of quarks). A summation over 
all quark flavors is understood. Using Eq. J2J we can write the above equation 
in the \i basis. Of interest to us here is the part 

L eff = 9 V2{(tan6 w (T 3 - Q)C« - T 3 C*)q LX i R qL 

+ tanOwCuQquXiLqii} (17) 
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Fig. 2. The LSP-quark interaction mediated by s-quark exchange. 



The above interaction is almost diagonal in the quark flavor. There exists, 
however, mixing between the s-quarks q~L and q R (of the same flavor) i.e. 



with 

tan.20u = 



q L = cosdqqi + sin9qq 2 , q R = -sinOq-qi + cos9qq 2 (18) 
m u (A + ficot/3) _ m d {A + fitanj3) 



1 •> ' m 2 dR -m 2 r +m%cos2p/2 



m u L ~ m I R + m 1cos2j3 /2 
Thus Eq. JT7|l becomes 

L e ff = gV2{[B L cos9 q q~LXiR ~ B R sin9 q q R xiL\qi 



(19) 



with 



+ [B L sin9 q q L xiR. + B R cos9 q q R xiL\ <?2j 

B L (q) = -iCf^an^ - \C^, q = u {charge 2/3) 
B L {q) = -\C$tan9 u + \Cjk, q = d (charge -1/3) 
B R (q) = \tan9 L0 C\\, q — u (charge 2/3) 
B R (q) = —\tan9 iA jC 11 : q — d (charge —1/3) 
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The effective four fermion interaction takes the form 

L eff = {9^) 2 {{B L cos9qq L xiR - B R sin9qq R xiL) 
-5—^ (B L cos8 q xiRqL - B R sin6qXiLqR) 

+ (B L sin6 q q L xiR + cos6gqRXlL)-n— 

Q - mq2 

(B L sin8qXiRqL + B R cos9qXiLqR.)} (20) 
The above effective interaction can be written as 

r _ T LL + RR , T LR / 01 \ 
L eff ~ L eff +L eff V 21 ) 

The first term involves quarks of the same chirality and is not much effected 
by the mixing (provided that it is small). The second term involves quarks of 
opposite chirality and is proportional to the s-quark mixing. 



The part L e ff 



LL+RR 



Employing a Fierz transformation L^ff can be cast in the more convenient 
form 



L^+ RR = (gV2) 2 2(~){\B L \ 2 

sin 2 £ 



J eff 



2' 



{-2 r -5 HLlxqLXlRl XlR 

q - m q 2 q z - rriq2 

. l9 , sin 2 9g cos 2 9a . x , ,„„. 

+ \B R \ 2 (- q - + - q —)qR.lxqRXiLl X XiL] (22) 

q - m S 2 q z - m q 2 

The factor of 2 comes from the Majorana nature of LSP and the (-1/2) comes 
from the Fierz transformation. Equation (12211 can be written more compactly 
as 

G F 

Leff = --j=2{qjx(PoR + 03RT3)(1 + lb)q 



with 



- qixiPoL + /3 3 lt 3 )(1 - 7sM(xi7 TXi} (23) 

Por = (gXl R + gxf J \C n tan6 w \ 2 

= {\xl R ~ \x\^) \C n tan6 w \ 2 (24) 
Pol = \\c*tane w + \c*\ 2 X l L + \\c*tan6 w - \c*\ 2 X \ 



'6 11 2 Zil ^ '6 11 2 

\~C*tan9 w + ~C*\ 2 xl L - \~C*tan9 w - i< 
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with 



X « L ~ « m 5 2 - g 2 + « - g 2 ' X «« ~ 9 m ~2 - g 2 + « 777,2 - g 2 1 j 
*i "12 **i ^2 

where Cg = cosOq, Sq = sinOq. The above parameters are functions of the 
four-momentum transfer which in our case is negligible. 
Eq H23|) it is often written as: 

L eff = -^2 [u lx (d°(u) + l5 d(u))u + d 7 x(d°(d) + l5 d(d))d] (xitVxi} 

(26) 

where 

d°(u) = p 0R + P 3R - Pol - Psl, d{u) = p 0R + p 3R + Pol + p 3L (27) 

d°(d) = fan - p 3R - Pol + P 3L ,d(u) = p 0R - P 3R + Pol - P 3 l (28) 

Proceeding as in sec. l3.ll we can obtain the effective Lagrangian at the nucleon 
level as 

L L ef L ; RR = -^(xxtVxi) Jx(q) (29) 
J X (q) = NlM(q) + ti(q)r 3 + f A {q)^ + f A (qh 5 T 3 )N (30) 

with 

fv = HPor - Pol), fv = 2gv(P 3 R ~ P 3 l) 
f° A = 2g A (p QR + Pol), f\ = 2g A (P 3R + P 3L ) (31) 

with g v = 1.0 and g A = 1.25. The quantity g A depends on the quark model 
for the nucleon. It can be anywhere between 0.12 and 1.00 (see below PPf . 

We should note that this interaction is more suppressed than the ordinary 
weak interaction by the fact that the masses of the s-quarks are usually larger 
than that of the gauge boson Z°. In the limit in which the LSP is a pure bino 
(Cn = 1, C21 = 0) we obtain 

Po R = tan 2 9 w {^xl R + lx} R ) , PtR = tan 2 9 w (J X 2 1J? - |x| R ) (32) 

Pol = (xl L + X\ L ) , Psl = — ^ (xl L - x\ ) (33) 

Assuming further that xa R — Xd R ~ Xu L — Xg L we obtain 

fHq)/fv(q) * +\ , fA(q)/f°A(q) * +^ (34) 

If, on the other hand, the LSP were the photino (Cn = cosBw,Ci\ = 
sin8yy,C 3 i = C41 = 0) and the s-quarks were degenerate there would be 
no coherent contribution (/° = if Pol = Por)- 

As we have mentioned in the previous section, to leading order, only the 
axial current contributes to the direct detection of the neutralino. 
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The part L^f f 
From Eq. ifflTjl we obtain 

Iff, = -(gV2) 2 sin26gB L ( q )B R (q)l[^ ^ ] (35) 

2 f- mg2 q z - mg2 

(QLXiRXiLqR + ^RXILXIRMl) 
Employing a Fierz transformation we can cast it in the form 

L eff = [(qqXiXi +9759X175X1 - (f^'ZXxio^Xi))] (36) 

where 

2 11 

(3{u) = -tan9 w C 11 {2sin20 il [-CR 11 tan9 w + -CR 2 i]A a (37) 

3 6 2 

/3(d) = sin26 s [^CR n tan6 w - ^C^}A^} (38) 

Where in the last expressions u indicates quarks with charge 2/3 and d quarks 
with charge -1/3. In all cases 

(m? -ml )ML 

— i/7 



and an analogous equation for Z\j. 

The appearance of scalar terms in s-quark exchange has been first 
noticed by Griest [3^1 ■ It has also been noticed there that one should consider 
explicitly the effects of quarks other than u and d m going from the quark 
to the nucleon level. We first notice that with the exception of t s-quark the 
qi — qn mixing small. Thus 

sin26,Au^ r 2m / A +^f )m \ , rtrMjM* ^M^tan(3)rn 2 w 
(mi L - q 2 ){mi R - q 2 ) " {mj - q 2 )(mj - q 2 ) 

(39) 

In going to the nucleon level and ignoring the negligible pseudoscalar and 
tensor components we only need modify the above expressions for all quarks, 
with the possible exception of the t quarks, by the substitution m q — > f q m^ 
(see sec. I4.1f) . For the t s-quark the mixing is complete, which implies that 
the amplitude is independent of the top quark mass. Hence in the case of the 
top quark we may not get an extra enhancement in going from the quark to 
the nucleon level. In any case this way we get 

L eff = ^[f s (q)NN + fh(q)Nr 3 N}xiXi (40) 

with 
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fs(q) 



U(3(u) + f d (3(d) 



(41) 



2 



q—s,c,b,t 



(42) 



2 



(see sec. 14. II for details). In the allowed SUSY parameter space considered in 
this work this contribution can be neglected in front of the Higgs exchange 
contribution. This happens because for quarks other than t the s-quark mixing 
is small. For the t-quark, as it has already been mentioned, we have large 
mixing, but we do not get the advantage of the mass enhancement. 

3.3 The Intermediate Higgs Contribution 

The coherent scattering can be mediated via the intermediate Higgs particles 
which survive as physical particles (see Fig. [3J . The relevant interaction can 



Fig. 3. The LSP-quark interaction mediated by Higgs exchange. 



arise out of the Higgs-Higgsino-gaugino interaction which takes the form 




+ H.C, 






(43) 
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Proceeding as above we can express W an B in terms of the appropriate 
eigenstates and retain the LSP to obtain 

- - tan9 w C^)C 3lX iBXiLHr) + B.C. (44) 



We can now proceed further and express the fields H\ , H® in terms of 
the physical fields h, H and A. The term which contains A will be neglected, 
since it yields only a pseudoscalar coupling which does not lead to coherence. 

Thus we can write 

£ eff = -^L X xN[fs(H) + fl(H)r 3 }N (45) 



fs ( H ) = o (ft* + 9d) + 9s + 9c + 9b + 9t (46) 



with 



where 

1 

2< 

fh{H) = \{9u-9d) (47) 

r ,,.cosa , „. sina-, m a . . , 

9 ai =[gi(h)—-+g 2 {H)—-]^, ai =u,c,t 48 
smp smfi tun 

9m = -Si" -s+92{H) , Kl = d,s,b (49) 

cosp cosp mjv 

c?i (fc) = A(C^tane w - C%)(C A1 cosa + Casino) m " mw ( 50 ) 

S 2 (#) = 4(C7f 1 tan0 w - CSfCusinct - C 3 icosa) WAfWw (51) 

m H — q A 

where mjy is the nucleon mass, and the parameters mh, niu and a depend 
on the SUSY parameter space (see Table 1). 



4 Going from the Quark to the Nucleon Level 

As we have already mentioned, one has to be a bit more careful in handling 
quarks other than u and d. 

4.1 The scalar interaction 

As we have seen the scalar couplings of the LSP to the quarks are proportional 
to their mass |19j . One encounters in the nucleon not only sea quarks (uu,dd 
and ss) but the heavier quarks also due to QCD effects J37]. This way one 
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obtains the scalar Higgs-nucleon coupling by using effective parameters f q 
defined as follows: 

(N\m q qq\N) = f q m N (52) 

where is the nucleon mass. The parameters /„, q = u,d, s can be obtained 
by chiral symmetry breaking terms in relation to phase shift and dispersion 
analysis. The isoscalar component can be obtained by considering the follow- 
ing quantities : 

1. The phenomenologically determined mass ratios: 



777 777 

— = 0.553 ±0.043 , — = 18.9 ±0.08 (53) 
m d m d 



2. The quantities : 

B u - B s 2B S 

z = m 1.49 , y = (54) 

Bd — B s Bd + B u 

One then finds that: 

B u 2z-{z-\)y B u 2 + (z-l)y 

W d = 2+(z-l)y f ° r pr0t ° nS ' -Bd = 2z~(z~ l)y ** nCUtr ° nS 

(55) 

with B q =< N\qq\N > 

3. The pion- nucleon sigma-term, ct^n'- 

this term is obtained from the isospin even 7T-N scattering amplitude with 
vanishing external momenta. It is defined by the scalar quark density 
operator averaged over the nucleon or equivalently by a^{t = 0), the 
scalar form factor of the nucleon at zero momentum transfer squared. The 
value of the sigma term is deduced from the analysis of two quantities: 
a-nNit — 2Af^) the scalar form factor at the Cheng-Dashen point, which is 
experimentally accessible, and the difference A a = cr 7r Ar(2M^) — ct w jv(0) = 
15.2 ± 0.4 MeV |2HI,[32I as induced by explicit chiral symmetry breaking. 
Experimentally, after efforts of many years, the value of the sigma-tern is 
still quite uncertain |39| . The canonical value of the ttN sigma term with 

777,, + 777rf , , , „ 

ovjv = {Bu + B d ) = (45± 8) MeV (56) 

is deduced from an earlier analyses with a T jv(t = 2M%) — 60 ± 8 [3"%| . 
During the last few years analyses of also more recent pion-nucleon scat- 
tering data lead to an increase in the value of scalar form factor at the 
Cheng-Dashen point a„ N (M%) with 88 ± 15 MeV gOj, 71 ± 9 MeV gT], 
79 ± 7 MeV [H] and (80 - 90) MeV gSj. Thus the recent analyses suggest 
a value for the pion-nucleon sigma term of about 

o-w = mu t md (B u + B d ) = (56 - 75) MeV (57) 
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4. Theoretical analysis of the a„N term: 

In the context of chiral perturbation one can show that: 

a nN = , cr = (35 ± 5) MeV (58) 

1 - V 

Eqs. I|56|) and (|57|l together with Eq. H58fl will provide the range of varia- 
tion in the parameter y. Taking: 

m u = 5.1 MeV , m d = 9.3 MeV (59) 

together with y will in turn provide by Ea. l55l the range of variation of the 
ratio B u /Bd- The uncertainties in Eqs. I|56I57I58[1 provide a wide range in 
which the parameter y can vary. In other words the experimental and the- 
oretical uncertainties permit us, we will exploit the possible consequences 
of variation in y to SUSY dark matter detection. For ovjv we choose 
45, 55, 65 and 75 MeV to reflect the range of values set by Eqs. and 
(15711 . Thus from Eq. (|58|l we extract the corresponding y parameters with 
0.22 ±0.11, 0.36 ± 0.09, 046 ± 0.08 and 0.53 ± 0.07, respectively. Then we 
will combine these values with Eq.E^to get the desired values of f q given 
below. 

From the above analysis we get in the case of the proton: 

yg^-^ti + ^-yV (60) 

0.756 m p 2 + (z — l)y 

fl = 0.553 f d [ 2Z ~l Z -!h , ff = H V (61) 

2 + (z — l)y 0.756 m p 2 

5 = 5^-^ (62) 

0.756 m p 2 

In the case of the neutron our expressions are analogous, the ratio B u /Bd 
getting the inverse value. 

For the heavy quarks, to leading order via quark loops and gluon exchange 
with the nucleon, one finds: 

fa = 2/27(1- £,/*) 

There is a correction to the above parameters coming from loops involving 
s-quarks [37] . The leading contribution can be absorbed into the definition, if 
the functions gi(h) and §2{H) as follows : 



gi{h )-, gim + l s{ 2<-^Ml^)} 
82(B) - gx(h)[l + |(2^ + ^?^|)] 



for Q = c and t For the b-quark we get: 



On The Direct Detection of Dark Matter 



15 



92 {H) -> 9l (h)[l + \(2<- - 

In addition to the above effects one has to consider QCD effects. These effects 
renormalize the contribution of the quark loops as follows |37| : 

fQCD{q) = I l+-y"(l„) 

with 

f3(a s ) = fj[l + f a s n} , lm {a s ) = 2^ 

Thus 

f Q CD{q) = l+ 1 i^ 

The QCD correction associated with the s-quark loops is: 

f Q CD(q) = l + f^ 

The above corrections depend on Q since a s must be evaluated at the scale 
of mQ. 

It convenient to introduce the factor /qc.d (9) / /qcd(<z) into the factors gi(h) 
and g2(H) and the factor of fQCD(q) into the the quantities /q. If, however, 
one restricts oneself to the large tanfl regime, the corrections due to the s- 
quark loops is independent of the parameters a and (3 and significant only for 
the t-quark. 

For a more detailed discussion we refer the reader to Refs. We 
thus obtain the results presented in Tabled 

We notice that there exist differences between the proton and neutron com- 
ponents. These, however, cannot be taken as the sole contribution to isovector 
contribution, since all quantities were derived with isoscalar operators. So the 
isovector contribution will be discussed elsewhere. Here we will limit ourselves 
to the isoscalar component f q = (fjj + /™)/2 

4.2 The axial current contribution 

The amplitudes a p = f\ + f\ and a n = f% — f\ are defined by gl| : 

a N = 2_j d q Aq N (63) 

q— u,d 7 s 

2s^Aq N = (N\q^ To q\N) (64) 

where s M is the nucleon spin and d q the relevant spin amplitudes at the quark 
level obtained in a given SUSY model. 

The isoscalar and the isovector axial current couplings at the nucleon level, 
/a 1 /A' are obtained from the corresponding ones given by the SUSY models 
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Table 1. The parameters f% and fg (upper part) as well as and /q (lower part) 
for the twelve cases discussed in the text. 



#rp rp rp rp rp rp rn rn rn rn rn rn 
Jd Ju Js Jc .Jb Jt Jd Ju Js Jc Jb It 

1 0.026 0.021 0.067 0.098 0.104 0.161 0.037 0.014 0.066 0.098 0.104 0.161 

2 0.027 0.020 0.133 0.087 0.092 0.144 0.037 0.015 0.133 0.086 0.092 0.143 

3 0.028 0.020 0.199 0.075 0.080 0.126 0.036 0.015 0.199 0.075 0.080 0.126 

4 0.033 0.025 0.199 0.078 0.083 0.132 0.044 0.018 0.199 0.077 0.083 0.122 

5 0.034 0.024 0.265 0.068 0.072 0.117 0.044 0.019 0.265 0.067 0.172 0.117 

6 0.031 0.025 0.332 0.057 0.061 0.106 0.043 0.017 0.332 0.057 0.062 0.102 

7 0.040 0.028 0.331 0.061 0.065 0.109 0.051 0.022 0.331 0.060 0.065 0.109 

8 0.041 0.028 0.400 0.051 0.055 0.095 0.051 0.023 0.400 0.050 0.055 0.095 

9 0.047 0.028 0.470 0.041 0.047 0.081 0.051 0.023 0.400 0.050 0.055 0.095 

10 0.047 0.027 0.462 0.045 0.050 0.090 0.050 0.023 0.470 0.040 0.044 0.060 

11 0.048 0.032 0.532 0.036 0.040 0.076 0.058 0.027 0.532 0.035 0.040 0.076 

12 0.049 0.032 0.603 0.026 0.030 0.063 0.057 0.027 0.603 0.026 0.030 0.063 



at the quark level, f%(q), /A(s)> via renormalization coefficients g\, g\, i.e. 
f% = S^/aC?)' /A = 9a /A (<?)• The renormalization coefficients are given terms 
of Aq defined above 02], via the relations 

g A = Au + Ad + As = 0.77 - 0.49 - 0.15 = 0.13 , g\ = Au - Ad = 1.26 

We see that, barring very unusual circumstances at the quark level, the 
isoscalar contribution is negligible. It is for this reason that one might prefer 
to work in the isospin basis. 

5 The nucleon cross sections 

With the above ingredients we are in a position to evaluate the nucleon cross 
sections. 

• The scalar cross section. As we have mentioned this is primarily due to 
the Higgs exchange. 

< x o = ^\fs + fs\ 2 , < x o = rolfs ~ fl? (65) 

with 



ctq = —{G F m p ) 2 = 0.77 x 10" 38 cm 2 = 0.77 x 10~ 2 pb (66) 



1 

2tt' 

Since, however, the process is dominated by quarks other than u and d, 
the isovector contribution is negligible. So we can talk about the nucleon 
cross section. 

The proton spin cross section is given by: 

=3a |/ A + /A| 2 =3a |a p | 2 (67) 
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6 The allowed SUSY Parameter Space 

It is clear from the above discussion that the nucleon cross section depends: 

• The the quark structure of the nucleon 

The allowed range of the parameters f q may induce variations in the nu- 
cleon cross section as large as an order of magnitude. 

• The parameters of supersymmetry. 

This is the most crucial input. One starts with a set of parameters at 
the GUT scale and predicts the low energy observable via the renormal- 
ization group equations (RGE). Conversely starting from the low energy 
phenomenology one can constrain the input parameters at the GUT scale. 

The parameter space is the most crucial. In SUSY models derived from min- 
imal SUGRA the allowed parameter space is characterized at the GUT scale 
by five parameters: 

• two universal mass parameters, one for the scalars, mo, and one for the 
fermions, my 2 . 

• tan(3. 

• The trilinear coupling Aq (or mP t ole ) and 

• The sign of fi in the Higgs self-coupling /iHi H2. 

The experimental constraints are 

1. The LSP relic abundance (including co- annihilations): 

0.09 < f2 LSP h 2 < 0.22 (previous), 0.09 < Q L sph 2 < 0.124 (WMAP) 

2. the b -> s 7 constraint (CLEO, BELLE) 

2 x 10~ 4 <BR< 4.1 x 10~ 4 

3. The Higgs mass: > 114.1 GeV. This applies on the standard model 
Higgs. So For SUSY one must correct for factor sin 2 (a — (3) where a is 
the Higgs mixing angle. So this imposes limits on tan (3 

4. Limits on g s — 2 (e~,e + ) experiments (E821 at BNL) 

a M = ($„ - 2)/2 = (33.7 ± 11.2) x 10- 10 

yields (2a level): 

11.3 x 10~ 10 < Sa^SUGRA) < 56.1 x 10~ 10 

5. The fermion masses: 

m t (pole) = 175 GeV , m b (m b ) = 4.25 GeV => 
m b (m z ) = 2.888 GeV, m T (M z ) = 1.7463 GeV 

We are not going to elaborate further on this interesting aspect, since it will 
be covered by another contribution to these proceedings by A. Lahanas. 
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7 Rates 



The differential non directional rate can be written as 

dRundir = . m da(u,v)\v\ (68) 

m x ArriN 

where A is the nuclear mass number, p(0) w 0.3GeV/cm 3 is the LSP density 
in our vicinity, m is the detector mass, m x is the LSP mass and da(u, v) is 
the differential cross section. 

The directional differential rate, i.e. that obtained, if nuclei recoiling in 
the direction e are observed, is given by |15| : 

dRdir — ^ —. lul'O.e OOv.e) — da(u, v S( — ^——j= — v.e) (69) 

m x Amjv 27r [i r v\J2 

where 0(x) is the Heaviside function. 

The differential cross section is given by: 

da(u, v) == 2 ^ U by)2 [{ZsF{uf + E spin F ll («)] (70) 

where u the energy transfer Q in dimcnsionlcss units given by 

u = -p— , Q = [m p Ab}- 2 = MeV (71) 

Qa 

with 6 is the nuclear (harmonic oscillator) size parameter. F(u) is the nu- 
clear form factor and Fn(u) is the spin response function associated with the 
isovector channel. 

The scalar cross section is given by: 



f 1 — i ^ r \ 2 ~s A 2 



1 , fh 2Z-A " 2 
1 + TTT 4 



7F 



(since the heavy quarks dominate the isovector contribution is negligible). 
a N x° * s ^ e LSP-nucleon scalar cross section. The spin Cross section is given 
by:' 

^spin — ( 7 r) O v o Cspinj Cspin To (^3) 

A*r(P) P ' X 3(1 + 4) 2 

S(«) « 5(0) = [An {Q)f + 24«o(0)« 1 (0) + J7i(0)) 2 ] (74) 
/a /a 

The couplings (/^) and the nuclear matrix elements J?i(0) (J?o(0)) associ- 
ated with the isovector (isoscalar) components are normalized so that, in the 
case of the proton at u = 0, they yield Cspin = !• 
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Table 2. The static spin matrix elements for various nuclei. For 3 He see Moulin, 
Mayet and Santos 48 , 49 . For the other light nuclei the calculations are from 
DIVARI gJJ. For 73 Ge and 127 I the results presented are from Ressel et al [2T] (*) 
and the Finish group et al |22| (**). For 207 Pb they were obtained by the Ioannina 
team (+). gU. 





3 He 


is p 


29 Si 23 Na 


73 Ge 127 P 


127p* 


207p b + 


O)(0) 


1.244 


1.616 


0.455 0.691 


1.075 1.815 


1.220 


0.552 




-1.527 


1.675 


-0.461 0.588 


-1.003 1.105 


1.230 


-0.480 


n p (o) 


-0.141 


1.646 


-0.003 0.640 


0.036 1.460 


1.225 


0.036 


o n \o) 


1.386 


-0.030 


0.459 0.051 


1.040 0.355 


-0.005 


0.516 


(J-th 




2.91 


-0.50 2.22 








f-lexp 




2.62 


-0.56 2.22 








Vth( s P in ) 




0.91 


0.99 0.57 









Mea;p 



With these definitions in the proton neutron representation we get: 



(spin = |s'(0) (75) 



S(0) = 



(— ft n (o)) 2 + 2-^„(o)fl p (o) + n 2 p (o) 

(Xp dp 



(76) 



where f2 p (0) and f2 n (0) are the proton and neutron components of the static 
spin nuclear matrix elements. In extracting limits on the nucleon cross sections 
from the data we will find it convenient to write: 



J? 2 (01 O 

spin * u "p\"/i J 'n , iS \2 IT7\ 

CT p f x o Capin = -^—IV^P + YfV^^ I ( 77 ) 



In Eq. 177|) 5 the relative phase between the two amplitudes a p and a n . The 
static spin matrix elements are obtained in the context of a given nuclear 
model. Some such matrix elements of interest to the planned experiments are 
given in table [21 The shown results are obtained from Divari 0^1, Ressel et al 
(*) HI], the Finish group (**) E2 and the Ioannina team (+) [35], 07]. 
The spin ME are defined as follows: 



fip(0) = \I'^~T~ x jj I°'*(p)I j J> ' ^n(°) = Y^y^^ J J\a z (n)\ J J y 

(78) 

where J is the total angular momentum of the nucleus and a z — 2S Z . The spin 
operator is defined by S z (p) = Yli=i i- e - a sum over a U protons in the 

nucleus, and S z (n) = Y^iLi S z {i), i.e. a sum over all neutrons. Furthermore 

ft o (0) = Q p (0) + a n (0) , fl x (0) = fip(0) - « n (0) (79) 
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8 Expressions for the Rates 

To obtain the total rates one must fold with LSP velocity distribution and 
integrate the above expressions over the energy transfer from Q m in determined 
by the detector energy cutoff to Qmax determined by the maximum LSP 
velocity (escape velocity, put in by hand in the Maxwellian distribution), 
i.e. v esc = 2.84 vo with vq the velocity of the sun around the center of the 
galaxy (229 Km/s). 

For a given velocity distribution f(u'), with respect to the center of the 
galaxy, one can find the velocity distribution in the Lab {(v,ve) by writing 
v = v + ve , ve=v$+ vi, with v i the Earth's velocity around the sun. 

It is convenient to choose a coordinate system so that x is radially out in 
the plane of the galaxy, z in the sun's direction of motion and y = x x z. 

Since the axis of the ecliptic lies very close to the x, y plane (uj — 186.3°) 
only the angle 7 = 29.8° (see Fig. becomes relevant. Thus the velocity of 
the earth around the sun is given by 

ve = vqz + t>i( sina x — cosa cos^y + cosa sirvy z ) (80) 

where a is phase of the earth's orbital motion. The LSP velocity distribution 




Fig. 4. The galactic plane is perpendicular to the paper containing the sun's velocity. 
The normal to the two planes form an angle 7 = 7r/2 — 7~7r/6. The modulation 
is affected by the projection of the Earth's velocity along the sun's velocity. Thus 
the velocity of the detector relative to the center of the galaxy is 220 + 15 = 235 
km/s around June 3nd (when the maximum of the event rate is expected) and 
220 — 15 = 205 km/s around December 3 (minimum of the event rate). 



i(v') is not known. Many velocity distributions are employed. In the present 
work we will adopt the standard practice and assume it to be Gaussian: 
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JV) 



-(V'/Vq) 2 



(81) 



Since Vi <^v we will ignore, for the moment, the motion of the Earth. 
Then the total (non directional) rate is given by 



R = P{0) m / 



R — Rt(a, Qmin) 
Pr ^2 fJ^KtS A2 



(82) 



m x o Am p v fJ. r (p) 



spin j- 



The SUSY parameters have been absorbed in R. The parameter t takes care 
of the nuclear form factor and the folding with LSP velocity distribution 
jSOl IS S3 E2 ■ It depends on Q m im i-e. the energy transfer cutoff imposed 
by the detector and a = [fij-bvoV^l^ 1 ■ 
In the present work we find it convenient to re-write it as: 



R = K 



where 



K 



P(0) 



Ccoh(A, /l r (A))(Jp X o + C sp i n (A, /l r (A))<Jp P _^ C, S pin 

p(0) m 



(83) 



100 GeV : 



vV) 



160 10" 4 (pb) 



-v 1 



0.3GeVcm~ s lKg 280/fcms- 1 

(84) 



and 



c coh (A, Hr(A)) = 



'--'Spin 



100 GeV 



lOOGeV 



fl r (A) 



-i 2 



(lr(p) 
Hr{A) 



Hr{p) 



A t coh (A) 

'spin (A) 



A 



(85) 



(86) 



The parameters c co h{A, /j, r (A)), c sp i n (A, n r (A)), which give the relative merit 
for the coherent and the spin contributions in the case of a nuclear target 
compared to those of the proton, have already been tabulated j^U for energy 
cutoff Qmin — 0, 10 keV. 

Via Eq. H83|) we can extract the nucleon cross section from the data (see 
below) . 

Neglecting the isoscalar contribution and using Q\ — 1.22 and Q\ = 2.8 for 
127 I and 19 F respectively the extracted nucleon cross sections satisfy: 



a p,x° 



C coh (A, fl r (A)) 
Cspin{A, H r (A)) 



xlO 4 (A = 127) , 



xlO 2 (A = 19) 



(87) 

It is for this reason that the limit on the spin nucleon cross section extracted 
from both targets is much poorer. 

The factors cl9 = c coh (19, /v(19)), sl9 = c spin (19, ^(19)), cl9 = c coh (73, ^(73)), 
s73 = c spm (73, Mr(73)) and cl27 = c co/l (127, Mr(127)), sl27 = c spi „(127, Mr(127)) 
for two values of Q m in and s3 = c sp , n (3, /x r (3)) for Q m in = can be found 
elsewhere |52| . 
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9 Bounds on the scalar proton cross section 

Before proceeding with the analysis of the spin contribution we would like to 
discuss the limits on the scalar proton cross section. In what follows we will 
employ for all targets H31-E1 the limit of CDMS II for the Ge target i.e. 
< 2.3 events for an exposure of 52.5 Kg-d with a threshold of 10 keV. This 
event rate is similar to that for other systems [KB] . The thus obtained limits 
are exhibited in Fig. [SJ 




o\03 








. 025 








0.02' 








.015 








.005 


80 


100 120 140 


160 180 200 




m.y — » GeV 



Fig. 5. The limits on the scalar proton cross section for A= 127 on the left and 
A= 73 on the right as functions of m x . The continuous (dashed) curves correspond 
to Qmin = (10) keV respectively. Note that the advantage of the larger nuclear 
mass number of the A= 127 system is counterbalanced by the favorable form factor 
dependence of the A= 73 system. 



10 Exclusion Plots in the a p , a n and cr p , cr n Planes 



From the data one can extract a restricted region in the a p , a n plane, which 
depends on the event rate and the LSP mass. Some such exclusion plots 
have already appeared |5b | -|57 | . One can plot the constraint imposed on the 
quantities \a p + jpa n \ and|y / o^' + 
limits via relations: 



a n e lS \ 2 derived from the experimental 



s2 n 



I V°p + jj-V°n) el I ^ cr bound (A) r(m x ,A), 



Abound (A) 



R 3 



10~ 5 pb 



K Ql cl™ n (A,Hr(A)) 



r(m x ,A) 



(A, Hr{A)) 



where 5 is the phase difference between the two amplitudes and clp® n (A, n r (Aj) 
is the value of c sp i n (A, fi r (A)) evaluated for the LSP mass of 100 GeV. Fur- 
thermore 



l a p + 7T" a «l ^ abound(A) [r(m x ,A)} 1/2 , a bound (A) = 

2 Irt 



{A) 

3(To 



1/2 



(89) 
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The constraints will be obtained using the functions cl°® n (A, fi r (A)), obtained 
without energy cut off , Q mi „ = 0, even though the experiments have energy 
cut offs greater than zero. Furthermore even though we know of no model 
such that e lS is complex, for completeness we will examine below this case 
as well. Such plots depend on the relative magnitude of the spin matrix el- 
ements. They will be given in units of the A-dependent quantity (Jbound(A) 
for the nucleon cross sections and the dimensionless quantity abound for the 
amplitudes respectively. Before we proceed further we should mention that, if 
both protons and neutrons contribute, the standard exclusion plot, must be 
replaced by a sequence of plots, one for each LSP mass or via three dimen- 
sional plots. We found it is adequate to provide one such plot for a standard 
LSP mass, e.g. 100 GeV, and zero energy threshold. The interested reader can 
find the scale for any other case in work already published JS2J ■ The situation 
is exhibited in Figs 16181 in the interesting case of the A=127 system using the 
nuclear matrix elements of Ressel et al given in Tabled For other targets we 
refer to the literature |52j . 

One can understand the asymmetry in the plot due to the fact that fl v is 
much larger than L?„. In other words if a p happens to be very small a large 
a n will be required to accommodate the data. In the example considered here, 
however, the extreme values differ only by 20% from the values on the axes, 
which arise, if one assumes that one mechanism at a time (proton or neutron) 
dominates. 




Fig. 6. The boundary in the a p ,a n plane extracted from the data for the target 
127 1 is shown assuming that the amplitudes are relatively real. The scale depends 
on the event rate and the LSP mass. Shown here is the scale for m x = 100 GeV. 
Note that the allowed region is confined when the amplitudes are of the same sign 
(left plot), but they are not confined when the amplitudes are of opposite sign. The 
allowed space now is i) The small triangle and ii) The space between the two parallel 
lines and on the right of the line that intercepts them. We also indicate by a dot 
the point a p — —a n favored by the spin structure of the nucleon. The nuclear ME 
employed were those of Ressel and Dean (see table 
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a p -> 5.0 x 1CT 3 pb 



Fig. 7. The same as in Fig. [5] for the a p ,a n plane. On the left the allowed region 
is that below the curve (the amplitudes are relatively real and have the same sign) . 
In the plots on the right the amplitudes are relatively real and of opposite sign. The 
allowed region is i) between the higher segment of the hyperbola and the straight 
line and ii) Between the straight line and the lower segment of the curve. The nuclear 
ME employed were those of Ressel and Dean (see table |5J 




b 0.2 o.4 o.g o.8 ia p _> 5.0 x 10~ 3 pb 

Fig. 8. The same as in Fig. [7| assuming that the amplitudes are not relatively real, 
but are characterized by a phase difference S. The allowed space is now confined. 
The results shown for the thin solid, thick solid and dashed curves correspond to 
5 = 7r/3, 7t/6 and ir/2 respectively . 



11 The modulation effect. 

As we have mentioned the expected event rate is so low that, even if one goes 
underground, the background is formidable. Especially since the signal coming 
from the detection of the energy energy of the recoiling nucleus has the same 
shape as that of the background. One, therefore, looks for specific signatures 
associated with the reaction. Since the event rate depends on the relative 
velocity between the LSP and the target, a periodic seasonal dependence is 
expected due to the motion of the Earth around the sun. What counts is the 
the is the projection of the velocity of the earth on the sun's velocity (see Fig. 

E}. 

If the effects of the motion of the Earth around the sun are included, the total 
non directional rate is given by 
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R — K [c coh (A, /i r (A))cr^ x o(l + h(a, Q min )cosa 



)] 



(90) 



and an analogous one for the spin contribution, h is the modulation amplitude, 
which is quite small, less than 2% and it depends on the velocity distribution, 
the nuclear form factor and, for a given target, on the LSP mass, a is the 
phase of the Earth, which is zero around June 2nd. In the case of the target 
127 I the modulation amplitude is shown in Fig.|5J We see that the modulation 



Fig. 9. The modulation amplitude h as a function of the LSP mass in the case of 
1271 for Qmin — on the left and Qmin = 10 keV on the right. We should mention 
that the average LSP energy for an LSP mass m x — 100 GeV is ~ 40 keV. For the 
definitions see text. 



amplitude is small, especially for Q m in = 0. Furthermore its sign is uncertain, 
since it depends on the LSP mass. The modulation amplitude increases as 
the threshold cut off energy increases, but, unfortunately, this occurs at the 
expense of the total number of counts. Furthermore many experimentalists 
worry that there are may be seasonal variations in the relevant backgrounds 
as well. 



12 Transitions to excited states 

As we have mentioned the average neutralino energy scales with its mass. It is 
~ 40 keV for m x = 100 GeV. Thus the neutralino energy is not high enough 
to excite the nucleus. In some rare cases involving odd mass nuclei there exist 
excited states at low energies, which can be populated in the LSP-nucleus 
collision due to the high velocity tail of the neutralino velocity distribution. 
From an experimental point of view this is very interesting |30| . since the 
signature of the 7— ray emission following the nuclear de-excitation is much 
easier than nuclear recoils. An interesting target is 127 I, which has an excited 
state at ~ 50 keV. It has recently been found |3J that the branching ratio to 
this excited state is appreciable from an experimental point of view. 




. 12 



m x -> (GeV) 
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rn x -» (GeV) 



Fig. 10. The ratio of the rate to the excited state divided by that of the ground 
state as a function of the LSP mass (in GeV) for 127 1. It was found that the static 
spin matrix element of the transition from the ground to the excited state is a factor 
of 1.9 larger than that involving the ground state. The spin response functions Fn(u) 
were assumed to be the same. On the left we show the results for Qmin = and 
on the right for Qmin = 10 KeV . In the last case, due to the detector energy cut, 
off the denominator (recoil rate) is reduced, while the numerator (the rate to the 
excited state) is not affected. 



13 The directional rates 

As we have already mentioned one may attempt to measure not only the en- 
ergy of the recoiling nucleus, but observe its direction of recoil. Admittedly 
such experiments are quite hard [221 , but they are expected to provide un- 
ambiguous signature against background rejection. Since the sun is moving 
around the galaxy in a directional experiment, i.e. one in which the direction 
of the recoiling nucleus is observed, one expects a strong correlation of the 
event rate with the motion of the sun In fact the directional rate can be 
written as: 

Rdir = ^—R t[l + h m cos(a - a m w)] (91) 

where h m is the modulation and a m is the "shift" in the phase of the Earth 
a, since now the maximum occurs at a = a m n. k/ (2tt) is the reduction factor 
of the unmodulated directional rate relative to the non-directional one. The 
parameters k , h m , a m depend on the direction of observation: 

e = (sin (9 cos ^ , sin (9 sin , cos 6*) 

The parameter Kt for a typical LSP mass 100 GeV is shown in Fig. ^2 as a 
function of the angle for the targets ^4=19 and A = 127. We see that the 
change of the rate as a function of the angle for the Maxwellian LSP velocity 
distribution is quite dramatic. This figure is important in the analysis of the 
angular correlations, since, among other things, there is always un uncertainty 
in the determination of the angle in a directional experiment. We prefer to 
use the parameters k and h m , since, being ratios, are expected to be less 
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Fig. 11. The quantity Kt as a function of the angle 0, the polar angle from the 
sun's direction of motion, for A — 19 on the left and A = 127 on the right. The 
results presented correspond to an LSP mass of 100 GeV. 



dependent on the parameters of the theory. We exhibit the dependence of the 
parameters t, h, k, h m , and a m , which are essentially independent of the LSP 
mass for target A = 19, in Table [3] (for the other light systems the results are 
almost identical). 

The asymmetry is quite large. For a Gaussian velocity distribution we find: 

As = „ o.97 

R(-z) + R(+z) 

In the other directions it depends on the phase of the Earth and is equal 
to almost twice the modulation. For a heavier nucleus the situation is a bit 
complicated. Now the parameters k and h m depend on the LSP mass |45| . 
It is clear that, if such experiments will ever be performed, such signatures 
cannot be mimicked by background events. 



14 Observation of electrons produced during the 
LSP-nucleus collisions 

Since the detection of recoiling nuclei is quite hard one may look for other 
events. One such possibility is the observation of ionization electrons produced 
directly during the LSP nuclear collisions Due to the properties of 

the bound electron wf, the event rate peaks at very low electron energies. 
One therefore must be able to achieve very low energy thresholds. In order to 
avoid uncertainties arising from the constraint SUSY parameter space we have 
opted to present the ratio of the event rate for producing electrons divided by 
the standard coherent recoil rate. This ratio is exhibited as a function of the 
electron threshold energy in Fig. El We see that for large atomic number Z 
and sufficiently low threshold energy this ratio may exceed unity. 

It has also been found that inner Is electrons can be ejected with a non 
negligible probability The produced electron holes can be filled via the 
Auger process or a sizable fraction can proceed via very hard (32 keV) X- 
ray emission. The detection of such X-rays, in or without coincidence with 
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Table 3. The parameters t, h, K,h m and a m for the isotropic Gaussian velocity 
distribution and Q m in = 0. The results presented are associated with the spin 
contribution, but those for the coherent mode are similar. The results shown are 
for the light systems. For intermediate and heavy nuclei there is a dependence on 
the LSP mass. +x is radially out of the galaxy (0 — n/2,<l> — 0), +z is in the 
sun's direction of motion (0 = 0) and +y is vertical to the plane of the galaxy 
(0 — tt/2,<P — it/2) so that (x,y,z) is right-handed. a m — 0,1/2,1,3/2 means 
that the maximum occurs on the 2nd of June, September, December and March 
respectively. 



type t h dir k h m en 



+z 0.0068 0.227 1 
dir +(-)x 0.080 0.272 3/2(1/2) 

+ (-)y 0.080 0.210 (1) 

-z 0.395 0.060 

all 1.00 

all 0.02 



nuclear recoils, will provide a signature very hard to miss, if SUSY allows for 
detectable recoil rates. 




-1 , 1 , 1 , 1 , o -I , 1 , i , ==f 

0.000 0.004 0.008 0.012 0.4 0.8 1.2 1.6 

T [keV] E th [keV] 



Fig. 12. On the left we show the differential rate for ionization electrons, divided by 
the total rate associated with the nuclear recoils, as a function of the electron energy 
T (in keV) for various atoms. On the right we show the total rate for producing 
electrons divided by the corresponding rate for nuclear recoil as a function of the 
threshold energy. The event rate is per atom, i.e. all electrons in the atom have been 
considered. The results exhibited were obtained for a typical LSP mass m x — 100 
GeV. 
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In this review we have dealt with various issues involving the direct detection 
of supersymmetric dark matter. The standard experiments employ various 
techniques of measuring the energy of the recoiling nuclei after their elastic 
scattering with the dark matter candidates. We have seen that the evaluation 
of the event rates involves a number of issues: 1) A supersymmetric model 
with a number of parameters, which at present can only be constrained from 
laboratory data at low energies as well as cosmological observations. 2) The 
dependence of the nucleon cross section on quarks other than u and d. 3) A 
proper nuclear model, which involves the nuclear form factor in the case of the 
the scalar interaction and the spin response function for the axial current. 4) 
Information about the density and the velocity distribution of the neutralino 
(halo model). 

Using the present experimental limits on the event rate and suitable inputs 
in 3)-4) we have derived constraints in the nucleon cross sections. Since the 
obtained event rates are extremely low, we have examined some additional 
signatures inherent in the neutralino nucleus interaction, such as the periodic 
behavior of the rates due to the motion of Earth (modulation effect). Since, 
unfortunately, this is characterized by a small amplitude, we were lead to 
examine the possibility of directional experiments. Tcsc, in addition to the 
recoil energy, will also attempt to measure the direction of the recoiling nuclei. 
The event rate in a given direction is ~ 6tt smaller than that of the standard 
experiments, but one maybe able to exploit two novel characteristic signatures: 
a) large asymmetries and b) interesting modulation patterns. 
Proceeding further we extended our study to include evaluation of the rates 
for other than recoil searches such as: i) Transitions to excited states and 
the observation of de-excitation 7 rays, ii) detection of the recoiling electrons 
produced during the neutralino-nucleus collision and iii) observation of hard 
X-rays, following the de-excitation of the ionized atom. 

With all the above signatures one hopes that, if the supersymmetric models 
do not conspire to lead to large suppression of the amplitudes, the direct 
direction of dark matter may soon follow. 

Acknowledgements 

This work was supported by European Union under the contract MRTN-CT- 
2004-503369 as well as the program PYTHAGORAS- 1. The latter is part of 
the Operational Program for Education and Initial Vocational Training of the 
Hellenic Ministry of Education under the 3rd Community Support Frame- 
work and the European Social Fund. The author is indebted to Professor Lef- 
teris Papantonopoulos for support and hospitality during the Aegean Summer 
School. 



30 J.D. Vergados 



References 

1. S. Hanary et at Astrophys. J. 545, L5 (2000); 
J.H.P Wu et al: Phys. Rev. Lett. 87, 251303 (2001); 
M.G. Santos et at Phys. Rev. Lett. 88, 241302 (2002). 

2. P.D. Mauskopf et at Astrophys. J. 536, L59 (2002); 
S. Mosi et at. Prog. Nuc.Part. Phys. 48, 243 (2002); 
S.B. Ruhl al, astro-ph/0212229 and references therein. 

3. N.W. Halverson et al: Astrophys. J. 568, 38 (2002) 

L.S. Sievers et at astro-ph/0205287 and references therein. 

4. G. Smoot et al (COBE Collaboration): Astrophys. J. 396, LI (1992). 

5. D. Spergel et al: Astrophys. J. Suppl. 148, 175 (2003). 

6. M. Tegmark et al: Phys.Rev. D 69, 103501 (2004). 

7. A. Jaffe et al: Phys. Rev. Lett. 86, 3475 (2001). 

8. G. Jungman, M. Kamionkowski, and K. Griest: Phys. Rep. 267, 195 (1996). 

9. D. Bennett et al.: Phys. Rev. Lett. 74, 2967 (1995). 

10. R. Bernabei et al: Phys. Lett. B 389, 757 (1996). 

11. R. Bernabei et al: Phys. Lett. B 424, 195 (1998). 

12. A. Benoit et al, [EDELWEISS collaboration]: Phys. Lett. B 545, 43 (2002); 
V. Sanglar,[EDELWEISS collaboration] arXiv:astro-ph/0306233 

D.S. Akerib et al,[CDM S Collaboration]: Phys. Rev D 68, 082002 (2003); 
arXiv:astro-ph/0405033 

13. G. Kane et al: Phys. Rev. D 49, 6173 (1994). 

14. A. Bottino et at, Phys. Lett B 402, 113 (1997). 

R. Arnowitt. and P. Nath, Ph ys. Rev. Lett. 74, 4952 (1995); Phys. Rev. D 54, 
2394 (1996); hep-ph/9902237 

V.A. Bednyakov, H.V. Klapdor-Kleingrothaus and S.G. Kovalenko, Phys. Lett. 
B 329, 5 (1994). 

15. J. Ellis, K. A. Olive, Y. Santoso, and V. C. Spanos: Phys.Rev. D 70, 055005 
(2004). 

16. M. W. Goodman and E. Witten: Phys. Rev. D 31, 3059 (1985). 

17. P. Ullio and M. Kamioknowski: JHEP 0103, 049 (2001). 

18. M.E.Gomez and J.D. Vergados, Phys. Lett. B 512 , 252 (2001); hep-ph/0012020 
M.E. Gomez, G. Lazarides and Pallis, C, Phys. Rev.D 61, 123512 (2000) and 
Phys. Lett. B 487, 313 (2000). 

19. M. Drees and N. N. Nojiri, Phys. Rev. D 48, 3843 (1993); Phys. Rev. D 47, 
4226 (1993). 

20. T.P. Cheng, Phys. Rev. D 38, 2869 (1988); H-Y. Cheng, Phys. Lett. B 219, 347 
(1989). 

21. M.T. Ressell et al, Phys. Rev. D 48, 5519 (1993); M.T. Ressell and D.J. Dean, 
Phys. Rev. C 56, 535 (1997). 

22. E. Homlund and M. Kortelainen and T.S. Kosmas and J. Suhonen and J. Toiva- 
nen, Phys. Lett B, 584,31 (2004); Phys. Atom. Nucl. 67, 1198 (2004). 

23. A. K. Drukier et al, Phys. Rev. D, 33, 3495 (1986); 
J.I. Collar et al., Phys. Lett B 275, 181 (1992). 

24. A. Green: Phys. Rev. D 66, 083003 (2002). 

25. P. Sikivie, I. Tkachev and Y. Wang, Phys. Rev. Let. 75, 2911 (1995); Phys. Rev. 
D 56, 1863 (1997) 

P. Sikivie, Phys. Let. b 432, 139 (1998); astro-ph/9810286 



On The Direct Detection of Dark Matter 



31 



26. D. Owen and J. D. Vergados: Astrophys. J. 589, 17 (2003); astro-ph/0203923 

27. J.D. Vergados and T.S. Kosmas, Physics of Atomic nuclei, Vol. 61, No 7, 1066 
(1998) (from Yadernaya Fisika, Vol. 61, No 7, 1166 (1998). 

28. K.N. Buckland, M.J. Lehner and G.E. Masek, in Proc. 3nd Int. Conf. on Dark 
Matter in Astro- and part. Phys. (Dark2000), Ed. H.V. Klapdor-Kleingrothaus, 
Springer Verlag (2000). 

29. The NAIAD ex periment B. Ahmed et al, Astropart. Phys. 19 (2003) 691; 
hep-ex/0301039 

B. Morgan, A. M. Green and N. J. C. Spooner, Phys. Rev. D 71 (2005) 103507; 
astro-ph/0408047 

30. H. Ejiri, K. Fushimi, and H. Ohsumi: Phys. Lett. B 317, 14 (1993). 

31. J. Vergados, PQuentin, and D. Strottman: IJMPE 14, 7 51 (2005). 

32. J. Vergados and H. Ejiri: Phys. Lett. B 606, 305 (2005); |hep^ph/040115l| 

33. C. C. Moustakidis, J. Vergados, and H. Ejiri: Nucl. Phys. B 727, 406 (2005). 

34. H. Ejiri and Ch.C. Moustakidis and J.D. Vergados, Dark matter search by ex- 
clusive studies of X-rays following WIMPs nuclear interactions, (to appear in 
Phys. Lett.); hep-ph/0507123 

35. J. D. Vergados: Part. Nucl. Lett. 106, 74 (2001); |hep^ph/0010151| 

36. K. Griest: Phys. Rev. Lett 61, 666 (1988). 

37. A. Djouadi and M. K. Drees, Phys. Lett. B 484, 183 (2000); S. Dawson, Nucl. 
Phys. B 359, 283 (1991); M. Spira it et al, Nucl. Phys. B453, 17 (1995). 

38. J. Gasser and H. Leutwyler and M.E. Sainio, Phys. Lett. B 253 (1991) 260; ibid 
B 253 (1991) 260. 

39. J. Gasser and H. Leutwyler: Phys. Rep. 87, 77 (1982). 

40. W. B. Kaufmann and G. E. Hite: Phys. Rev. C 60, 055294 (1999). 

41. M. G. Olsson: Phys. Lett. B 482, 50 (2000); |arXiT!hep -ph/0001203 

42. M. M. Pavan, I. I. Strakovsky, R. L. Workman, and R. A. Arndt: PiN Newslett. 
16, 110 (2002) ; |arXrv:hep-ph /01 1 1066 

43. M. G. Olsson and W. B. Kaufmann: PiN Newslett. 16, 382 (2002); 
ai'Xiv:hcp-ph/0111066 

44. The Strange Spin of the Nucleon, J. Ellis and M. Karliner, hep-ph/9501280 

45. J. D. Vergados: J.Phys. G 30, 1127 (2004). 

46. P. C. Divari, T. Kosmas, J. D. Vergados, and L. Skouras: Phys. Rev. C 61, 
044612 (2000). 

47. T. Kosmas and J. Vergados: Phys. Rev. D 55, 1752 (1997). 

48. E. Moulin, F. Mayet, and D. Santos: Phys. Lett. B 614, 143 (2005). 

49. D. Santos et al, The MIMAC-He3 Collaboration, A New 3 He Detector for non 
Baryonic Dark Matter Search, Invited talk in idm2004 (to appear in the pro- 
ceedings) . 

50. J. D. Vergados: Phys. Rev. D 62, 023519 (2000). 

51. J. D. Vergados: Phys. Rev. D 63, 06351 (2001). 

52. J.D. Vergados, Direct SUSY Dark Matter Detection- Constraints on the Spin 
Cross Section, hep-ph/0512305 

53. P. Belli, R. Ceru lli, N. Fornego, and S. Scopel: Phys. Rev. D 66, 043503 (2002); 
hep-ph/0203242 

54. M. M. Pavan, R.A. Arndt, I.I. Strakovsky and R.L. Workman, hep=ph/0111066. 

55. D. S. Akerib et al (CDMS Collaboration): Phys.Rev.Lett. 93, 211301 (2004). 

56. C. Savage, P. Gondolo, and K. Freese: Phys. Rev. D 70, 123513 (2004). 

57. F. Giuliani and T. Girard: Phys.Lett. B 588, 151 (2004). 



